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Abstract 

The Borromean nucleus ^^Ne {}^0+p + p) is investigated by using the hyperspheric 
adiabatic expansion for a a three-body system. The measured size of ^^O and the 
low-lying resonances of ^^F i^^O+p) are first used as constraints to determine both 
central and spin-dependent two-body interactions. Then, the ground state struc- 
ture of ^^Ne is found to be an almost equal mixture of s"^ and d? proton-^^0 relative 
states, the two lowest excited states have about 80% of sd-mixed components, and 
for the next two excited three-body states the proton-^^0 relative s-states do not 
contribute. The spatial extension is as in ordinary nuclei. The widths of the reso- 
nances are estimated by the WKB transmission through the adiabatic potentials 
and found in agreement with the established experimental limits. We compare with 
experimental information and previous works. 



PACS: 21.45.+V, 27.20. +n, 21.10.Dr 



1 Introduction 



Halo nuclei are weakly bound and spatially extended systems, and they are 
expected to appear along the driplines, where nucleon single particle s or 
p-states occur with sufficiently small separation energy [1-4]. Few-body tech- 
niques have proved to be very successful to describe them [5]. The core degrees 
of freedom decouple from the ones of the halo nucleons, and a cluster descrip- 
tion of the system becomes then appropriate [6] . 
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Among the halo structures, the three-body Borromean bound systems (where 
the two-body subsystems arc all unbound) are specially challenging [5,7,8]. 
For nuclei the most likely three-body halo candidates are Borromean systems 
with small binding energy. The reason is that three particles with attractive 
short-range interactions are more bound than only two of them. This is not 
necessarily true when also a sufficiently strong repulsive Coulomb interaction 
or Pauli repulsion is present. Furthermore the core should contain an even 
number of the nucleons surrounding it, because otherwise one valence nucleon 
would be pulled into the core leaving at most one with the larger spatial 
extension of a halo. 

Along the neutron dripline the most studied three-body Borromean halos are 
^He and ^^Li. Heavier neutron dripline systems present very interesting chal- 
lenges, but the interaction ingredients are often not well established or the 
simple three-body structures are not directly applicable. Examples could be 
^^Be, ^^B, ^^C, ^^F and ^^Ne. Since the Coulomb interaction works against halo 
formation, on the proton dripline the appearance of halos is less likely, and 
the candidates necessarily require a relatively light core. In fact, heavier halos 
than Ne arc not expected [9]. Still the degrees of freedom describing the core 
and the surrounding protons may decouple to some extent and a few-body 
treatment might be justified. 

The interest in genuine two-proton decays have increased recently especially 
due to the improved experimental techniques [10-12]. This fact has motivated 
different investigations on light proton rich nuclei. An especially interesting 
case is the lightest Borromean dripline nucleus ^''Ne (^^0-|-p+p), whose ex- 
cited resonance states are promising candidates to decay through direct two- 
proton emission [13,14]. However, the structure of ^^Ne is very controversial. 
The available publications predict very different content of the s and ci-waves 
[15-18]. Inconsistencies prevail even in a very recent and detailed three-body 
calculation of the structure of ^^Ne [19]. Nevertheless these models are used 
to predict new features of the delicate Thomas-Ehrman shifts between states 
in i^Ne and ^^N [20]. 

These inconsistencies may be related to the fact that ^^O has a non- vanishing 
spin which is difficult to handle consistently and therefore sometimes simply 
ignored. Unfortunately the finite core spin is a crucial ingredient in the struc- 
ture of the excited resonance states, at least if they have to be related to 
the measured resonances of ^^F nucleus. We therefore decided to investigate 
^^Ne with an established method which is well tested on the neutron dripline 
[21]. This method employs hyperspheric adiabatic expansion which is espe- 
cially suited to describe large and weakly bound systems like halo nuclei [22]. 
Furthermore, simultaneous use of the complex scaling transformation converts 
the method into an efficient tool to compute three-body resonances [23,24]. 
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The inconsistencies might perhaps be related to simphfied treatments of ^ Ne 
as a three-body system of a core and two protons. However, ah the lowest 
excited states in the ^^0-core occur with positive parities between 5.2 MeV 
and 9 MeV with one exception of 3/2~ at about 6.2 MeV. The ground, first 
and second excited states of -"^^Ne occur with negative parities at -0.94 MeV, 
0.34 MeV and 0.82 MeV from the two-proton threshold. These energies are 
small compared to even the lowest-lying core-excited states, which further- 
more would have to be combined with proton valence states of negative parity 
only present in the next shell. Thus two suppressing effects must occur simul- 
taneously to contribute to these states in ^'''Ne. In fact these arguments also 
apply to the third and fourth excited state in ^^Ne although now the energies 
are only a factor of two smaller than the core-excited states. 

The 3/2~ state at 6.2 MeV is then the most likely contributor, but its energy 
is still relatively high. We can not exclude small contributions form this core 
excitation but it seems unlikely that the ground state is affected, since then the 
two protons should pay a prize of combining into either one d-wave or two p- 
waves. Again this is double suppression and rather unlikely. Still contributions 
to the excited states of ^^Ne are not as strongly excluded. In any case such a 
core-excitation cannot be the origin of the inconsistencies. 

Insignificant contributions from core- excitations is equivalent to the assump- 
tion of a structureless core. The probability distributions of core and valence 
particles may still be overlapping provided the corresponding degrees of free- 
dom are dynamically decoupled. The restriction to the smaller Hilbert space 
in a few-body treatment necessitates a matching renormalization of the in- 
teraction achieved by phenomenological constraints on the effective potential. 
The related difficulty of antisymmetry between core and valence nucleons is 
approximately accounted for by use of phase equivalent potentials [25]. 

In the present work we investigate in detail the characteristics of the low- 
lying states in ^'^Ne in a three-body model. We arrive at a consistent picture 
describing simultaneously the two-body properties of the internal subsystems 
and the global three-body properties. The contributions to the ^^Ne states 
from the different partial waves are carefully analyzed. We do not encounter 
any sign of a break between properties of ground and excited states as a signal 
of substantial contribution of core-excitation from the possible 3/2~ state. 
The convergence of the model and possible uncertainties from the details of 
the two-body interactions are also investigated. In this attempt to arrive at a 
consistent picture we shall confront our calculations with results from previous 
work and if possible determine the origin of discrepancies. We already in [26] 
discussed a specific application to the Thomas-Ehrman shifts. 

The paper is organized as follows: In section 2 we describe very briefiy the 
method used to compute three-body bound state wave functions and three- 
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body resonances as well as the proton-proton and proton-^^0 interactions. 
With the three-body model and the input properly defined we then in section 
3 discuss the structure of the low-lying states of ^''Ne. The accuracy of the 
results and the relation to previous works are discussed in sections 4 and 5. 
We close the paper with a summary and conclusions. 



2 Bcisic ingredients 



An accurate description of a three-body system requires a suitable method 
to solve the Faddccv (or Schrodingcr) equations and the interactions between 
the constituent particles. We first very briefiy sketch the method and then we 
give a few details about the most important of the two-body interactions 



2.1 The three-body method 



We use the hyperspheric adiabatic expansion method as described in details 
for instance in [7,21]. The three-body wave function is a sum of three Faddeev 
components ip'^''\xi,y^) (i=l,2,3), each of them expressed in one of the three 
possible sets of Jacobi coordinates {ajj, t/j}. Each component is then expanded 
in terms of a complete set of angular functions {0^*^} 

P n 

where p — \/x^ -\- y^, ctj — arctan(xi/7/i), VLxi-i and ^ly^ are the hyperspheric 
coordinates. Writing the Faddeev equations in terms of these coordinates, and 
inserting the expansions (1), the Faddeev equations can be separated into 
angular and radial parts: 



n 



(2) 



^2 



dp' 



+ j^ivM -E) + j, [Up) + -j) 



dp 



-Qr 



Up) 

fn'{p)^0 



(3) 



where Vjk is the two-body interaction between particles j and k, is an 
angular operator [21] and m is the normalization mass. In Eq.(3) E is the 
three-body energy, V36 is a three-body potential used for fine-tuning and the 
functions P„„/ and Qnn' are given for instance in [21]. 
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It is important to note that the set of angular functions used in the expan- 
sion (1) are precisely the eigenfunctions of the angular part of the Faddeev 
equations. Their corresponding eigenvalues are denoted by A„(p), and enter in 
the radial equations (3) as an essential ingredient in the diagonal part of the 
effective radial potentials: 

T. / X Kip) + 15/4 , , 

VeS(p) = ^ p2 + Vsbip) (4) 

The complex scaling method [27] amounts to rotation of the hyperradius (p — > 
pe*^), while the hyperangles remain unchanged [24]. As soon as the rotation 
angle 9 is larger than the argument of any given resonance, then this resonance 
together with the bound states is obtained as an exponentially decreasing 
solution to the coupled set of radial equations (3). The problem of matching 
with the non trivial Coulomb asymptotics is thus avoided, since the resonance 
corresponds to a solution with "bound state" asymptotics. It was shown in 
[24] that the complex scahng method allows rehable calculations of resonance 
energies, also for three-body systems with the long-range Coulomb interaction. 

2.2 Two-body potentials 

The two-body potentials must describe the properties of the two-body sub- 
systems. However, unreflected reproduction of the low-lying spectra of the 
two-body systems is not sufficient. The subsequent use in the context of a spe- 
cific Hilbert space could be crucial. This becomes particularly evident when a 
number of degrees of freedom are frozen in the three-body computation while 
simultaneously also necessary to conserve symmetries of the effective Hamil- 
tonian. An appropriate example is the spin-dependence of the assumed form 
of the effective two-body potentials. A bad choice can lead to catastrophic 
results for the three-body system [28] . 

For the proton-proton interaction we use the parametrization given in [26,29]. 
Other choices have been tried as well but the three-body results do not change 
significantly. 

For the proton-^^0 interaction spin-dependent operators are indispensable to 
reproduce the properties of the low-lying ^^F resonances. We parametrize by 
use of operators that conserve the proton total angular momentum jp—i+Sp 
and the total two-body angular momentum j—jp+Sc, where i is the relative 
orbital angular momentum, and Sc and Sp are the core (^^O) and proton 
spins, respectively. These operators permit a clear energy separation of the 
usual mean-field spin-orbit partners ii+1/2 and it-i/2- It is then possible to 
use a proton-core interaction such that the low-lying states have well defined 
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Table 1 



Radius (6^) and strengths of the central (Si^^), spin-spin (S'ssO) spin-orbit 
{Sio) potentials in Eqs.(5) and (6), that correspond to the Woods-Saxon (WS) and 
Gaussian (G) proton-core interactions, respectively. In (5) the difFuseness a is 0.65 
fm in all the cases. 



£ 


be (fm) 


Si^^ (MeV) 


Si'^ (MeV) 


Si^j (MeV-fm2, MeV) 




WS G 


WS G 


WS G 


WS G 





3.00 2.60 


-53.91 -76.56 


0.92 1.34 




1 


2.70 2.80 


-19.99 -21.42 


0.69 1.04 


-25.0 -8.2 


2 


2.85 2.80 


-58.45 -75.89 


0.24 0.32 


-25.0 -8.2 



ij^ quantum numbers, like the ^5/2 states in F. This is especially important 
when one of these states is forbidden by the Pauli principle, like for instance 
the P3/2 waves in ^°Li [28]. 



In [26] we introduced an ^-dependent proton-core interaction of the form: 



core(^) (^) 

5f /f (r) + S^f^{r)s, . - Si:^l^^ni\r)£ ■ s, + ^Erf(r/6,) , 

where the shapes of the central (fP), spin-spin (/jf-*) and spin-orbit (f^^) ra- 
dial potentials are chosen to be Woods-Saxon functions, 1 / ( H-exp {{r — be)/a)), 
with the same diffuseness a in all terms. The proton number of the core is 
Zc- The error function Erf describes the proton-core Coulomb interaction of a 
gaussian core-charge distribution. The range parameter 6c=2.16 fm reproduces 
the experimental rms charge radius in ^^O of 2.65 fm. 

In [26] the parameters of the Woods-Saxon radial potentials were adjusted 
to reproduce the experimental energies of the four lowest states in ^^F, and 
simultaneously, after switching off the Coulomb interaction, the energies of 
the corresponding states in the mirror nucleus ^^N. The resulting parameters 
are given in table 1. 

In order to investigate the dependence of the results on the two-body interac- 
tions we have constructed an additional proton-core potential with gaussian 
form factors. Assuming that the (possibly ^-dependent) range of the gaus- 
sians is the same for central, spin-spin, and spin-orbit terms we can write the 
two-body interaction as: 

V^^Ur) - e-(^/^^)^ (5f + Si^Js, . 3^ + Si^Ji • . (6) 



The s and d-wave parameters are adjusted to reproduce the experimental 
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Table 2 

The two-body ^^N and ^^F spectra obtained with the gaussian potential specified in 
table 1. For unbound states we give the corresponding energies and widths {Er, T). 
The experimental error bars are not specified when they are smaller than the last 
digit. For unbound states the energies are decay energies above threshold. 





16 N 


Exp. [30] 


16p 


Exp. [30] 


0" 


-2.37 


-2.371 


(0.53,2 • 10-3) 


(0.535, 0.040 ±0.020) 


1- 


-2.09 


-2.094 


(0.70,~ 10-5) 


(0.728 ± 0.006, < 0.040) 


2- 


-2.49 


-2.491 


(0.95,0.01) 


(0.959 ± 0.005, 0.040 ± 0.030) 


3- 


-2.19 


-2.193 


(1.21,0.01) 


(1.256 ±0.004, < 0.015) 


1+ 


(0.86,1.22) 


(0.86,0.02) 


(3.71,3.80) 


(4.29±0.01, < 0.040) 


2+ 


(1.03,1.86) 


(1.03,<0.01) 


(3.89,4.31) 


(4.41±0.01, < 0.020) 


1- 


(2.33,0.96) 


(2.27±0.05,0.25±0.05) 


(5.28,2.27) 


(5.81±0.01,— ) 


2" 


(2.44,1.07) 


(2.56,0.02±0.01) 


(5.38,2.41) 





energies of the 0~, 1^, 2^, and 3^ states in ^^N and ^^F, see table 1. The 
accuracy is rather good and very similar to the Woods-Saxon case [28] as seen 
in table 2 where both computed and measured energies are given. 

The experimental energies of the second l~/2~ doublet is used to adjust the 
value of the d-wave spin-orbit parameter. These two states must come from 
the coupling of a d3/2-wave and the 1/2" state of the core. In the same way, 
the measured l~^/2'^ doublet, coming from the coupling of a p3/2-wave and the 
1/2^ state of the core, is used to determine the p-wave interaction (also given 
in table 1). As seen in the lower part of table 2, the computed resonance ener- 
gies agree reasonably well with measurements, but the widths are very much 
larger than the experimental values. The reason is simply that the computed 
resonances appear above the Coulomb and centrifugal barriers. This in turn 
is the result of the gaussian form of the attractive central potential combined 
with the necessary range of about 3 fm. Then the attraction at the surface 
pushes the barriers outwards and simultaneously reduces their heights. 

It is possible to construct other two-body miclcon-corc potentials reproducing 
as well the widths of the high-lying 1^/2^ and 1+/2+ doublets. This would 
require a more complicated form than Eq. (6) and more than the four parame- 
ters used for each partial wave. With the present assumptions we reproduce 6 
experimental values for s-waves and 10 for d-waves. We prefer the simphcity 
of Eq.(6), because components related to the relatively high-lying ^3/2 state as 
well as the p-waves in general, as we shall see later, never contribute by more 
than 4% of the total norm of the wave function. Another potential would es- 
sentially not change the results as long as symmetry properties and resonance 
energies are maintained, see the discussion of accuracy in section 4. 

The Woods-Saxon and Gaussian s-wave potentials in table 1 have a deeply 
bound ^^F-state at —26.2 MeV and —30.5 MeV, respectively. These states 
correspond to the Si/2 proton states occupied in the ^^O core. They are then 
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forbidden by the Pauli principle, and should be excluded from the calcTilation. 
This is implemented as in rcfs. [29,25] by use of the phase equivalent potential 
which has exactly the same phase shifts as the initial two-body interaction 
for all energies, but the Pauli forbidden bound state is removed from the two- 
body spectrum. We then use the phase equivalent potential of the central 
part of the s-wave potentials in table 1. Thus the s-states actively entering 
the three-body calculations are the second states of the initial potentials. 

The lowest p shell is also fully occupied in ^^O. We should then apply the 
same treatment as for s-waves to the p-wave proton-core interaction, using a 
potential with deeply bound states that are afterwards removed by the cor- 
responding phase equivalent potentials. However, since the p waves basically 
have no effects in the three-body calculation we use for simplicity the shallow 
l—l potential without bound states as given in table 1. This saves a substantial 
amount of computing time without loss of accuracy in the calculations. 

Partial waves with ^ > 2 are not crucial for the low-lying states. They should 
contribute on roughly the same level as the positive parity core-excitations 
corresponding to energy increases from the p to the sd-shell, and perhaps less 
than the negative parity excitation of ^^O increasing the energy only by the 
p-shcll spin-orbit splitting. Quantitative estimates can be obtained with po- 
tentials systematically extrapolated from the £ = 0, 1,2 parameters in table 
1. This can then be inferred to indicate that core-excitations are insignificant. 
In any case the use of effective potentials reproducing the pertinent measured 
quantities has to be consistent with the choice of model space. Then a phe- 
nomenologically adjusted two-body model without core- excitations for ^^F is 
expected to account correctly for most of the energy in the ^^Ne spectrum. 



3 Structure of the ^^Ne states 

In this section we discuss the results obtained for the ground and excited 
states of ^'''Ne. We start with the results obtained with the Woods-Saxon (WS) 
proton-core potential in Eq.(5). 

Together with the two-body potentials we must specify the components in- 
cluded in the calculations. We give them in subsection 3.1. The output are then 
the effective potentials entering in the radial part of the Faddeev equations, 
and the energy spectrum obtained when solving the radial part. The poten- 
tials are briefly described in subsection 3.2, and in 3.3 the properties of the 
computed -"^^Ne states are analyzed in detail. We close the section by discussing 
the spatial distribution of the states and provide a novel WKB estimate of the 
widths of the resonances in subsections 3.4 and 3.5, respectively. 
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Table 3 

Components included for the l/2~ state of ^''Ne. The upper part refers to the 
Jacobi set where x connects the two protons. In the lower part x connects the 
core and one of the protons. The sixth column indicates the maximum value of the 
hypermomentum K. The last column gives the contribution of the component to 
the total norm of the wave function. Only those components contributing more than 
1% are given. In each Jacobi set these probabilities should then add up to 100%. 
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44 
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1/2 


42 
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1/2 
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3. 1 Components 



We first compute the angular eigenvectors and eigenvalues from Eq.(2). This is 
done by expanding the eigenvectors in the basis {y^^i^^Li'^h ^yi)®Xs.^sy,s}, 
where y^^^y ^ are the hyperspheric harmonics and x is the spin function [21]. 
The quantum number £^ is the relative orbital angular momentum between 
particles j and A;, and Sx is the coupling of their two spins, is the relative 
orbital angular momentum of particle i and the center of mass of the jk two- 
body system, and Sy is the spin of particle i. The coupling of Ix and ly, and of 
Sx and Sy^ are L and S", respectively, which in turn couple to the total angular 
momentum J of the system. Finally, K=2n+lx+(!-yi where n is the usual nodal 
quantum number. 

The lowest valence space for the two protons in ^^Ne is the sci-shell. Therefore 
the expected dominating components are 1^—^,2, £y—0,2. By use of these 
components and the core-spin and parity of l/2~, we can construct three-body 
states with total angular momentum and parity J'^=l/2~, 3/2^, 5/2^, 7/2~, 
and 9/2~. Nevertheless, it is only after a full calculation that it is possible to 
determine precisely which components are needed to reach a given accuracy. 
For the states mentioned above the components are specified in tables 3 to 7. 
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Table 4 



Table 5 



*LS table 3 for the 3/2" 


state of I'^Ne. 


As table 3 for the 5/2" 


state of I'^Ne. 






L 




s 










£y 


L 


Sec 


s 




Vy(i^Ne) 


1 


1 


1 


1 


1/2 


62 


5.7 




1 


1 


1 


1 


3/2 


62 


6.8 


1 


1 


1 


1 


3/2 


42 


1.2 




1 


1 


2 


1 


1/2 


52 


4.2 
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1 


2 


1 


1/2 


42 
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1 


3/2 


62 


7.3 


1 


1 


2 


1 


3/2 


42 


5.0 
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2 


2 





1/2 


84 


9.1 
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1/2 


54 
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1 


1/2 


32 


2.1 


2 


2 


1 


1 


1/2 


54 


1.0 




2 


2 


1 


1 


3/2 


64 


5.8 


2 


2 


1 


1 


3/2 


34 


0.7 




2 


2 


2 





1/2 


64 


3.6 


2 


2 


2 





1/2 


54 


3.3 




2 


2 


2 


1 


1/2 


74 


10.6 


2 


2 


2 


1 


1/2 


84 


10.1 




2 





2 





1/2 


42 


3.1 





2 


2 





1/2 


42 


1.4 
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2 
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1/2 


82 


29.1 
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1/2 


100 


32.9 
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3/2 


62 
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1 


3/2 


54 


4.9 
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2 





1/2 


82 


16.6 


2 





2 





1/2 


94 


21.3 
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2 


1 


1/2 


82 


14.2 


2 





2 


1 


1/2 


74 


14.0 







2 


2 


1 


3/2 


62 


6.6 


2 





2 


1 


3/2 


44 


4.9 

















Table 6 

As table 3 for the 7/2" state of ^'^Ne. 



Table 7 





iy 


L 


Sx 


s 


Kmax 




2 


2 


4 





1/2 


44 


9.8 





4 


4 





1/2 


54 


19.4 


4 





4 





1/2 


54 


14.6 


1 


3 


3 


1 


1/2 


54 


20.2 


1 


3 


3 


1 


3/2 


44 


7.0 


3 


1 


3 


1 


1/2 


54 


21.5 


3 


1 


3 


1 


3/2 


44 


7.5 


2 


2 


3 





1/2 


124 


31.5 


2 


2 


3 


1 


1/2 


94 


10.6 


2 


2 


3 


1 


3/2 


94 


14.9 


2 


2 


4 





1/2 


64 


11.1 


2 


2 


4 


1 


1/2 


84 


31.4 



Is table 3 for the 9/2" 


state of i^Ne. 




£y 


L 




s 


^max 


VK("Ne) 


2 


2 


4 





1/2 


4A 


9.6 





4 


4 





1/2 


64 


18.9 


4 





4 





1/2 


64 


14.2 


3 


1 


3 


1 


3/2 


84 


29.3 


1 


3 


3 


1 


3/2 


84 


27.3 


2 


2 


3 


1 


3/2 


104 


56.0 


2 


2 


4 





1/2 


64 


11.2 


2 


2 


4 


1 


1/2 


84 


30.4 


2 


2 


4 


1 


3/2 


64 


1.9 



3.2 Effective potentials 



In Fig. la we show the deepest angular eigenvalue Ai obtained from equation 
(2) for each of the 1/2", 3/2", 5/2", 7/2", and 9/2" states in ^^Ne. Their main 
characteristics are the hyperspherical values K{K + 4) at p=0, a minimum due 
to the attractive short-range interaction and a linear increase at asymptotically 
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Fig. 1. a) The lowest angular eigenvalues Xi{p) for the l/2~, 3/2~, 5/2~, 7/2~ 
and 9/2~ states of ^''Nc as fvmction of p. b) Outer part: Three different effective 
three-body gaussian potentials that reproduce the experimental two-proton separa- 
tion energy of the l/2~ state in ^^Ne. Their parameters are given in table 8. The 
curve diverging at the origin is the corresponding effective potentials in Eq.(4). In- 
ner part: Comparison of the deepest effective potential (Eq. (4)) for the different 
J~ states of ^^Ne. The depth of the effective potentials decreases when J increases. 

large distances due to the repulsive Coulomb interaction [31]. The A-function 
starting at {K=0) is not present due to the use of the s-wave phase equivalent 
potential removing the deeply bound s-state to account for the Pauli principle 
[25]. Details about the higher A-functions are shown in [26]. 

From the figure we see that the most attractive A function in ^^Ne corresponds 
to the l/2~ state. This agrees with the fact that the l/2~ state is the only 
bound state in ^'''Ne with a measured two-proton separation energy of —944 
keV. Nevertheless, after the computation, the resulting ^'''Ne ground state bind- 
ing energy is —0.79 McV, i.e. 150 keV less bound than the experimental value. 
This is actually expected, since three-body calculations using pure two-body 
interactions typically underbind the system. This problem is solved by inclu- 
sion of the weak effective three-body potential V36 in Eq.(3), that accounts for 
three-body polarization effects arising when the three particles all are close to 
each other. Therefore the three-body potential has to be of short range, and 
its precise shape is not very relevant, since the overall three-body structure 
is unchanged. This construction ensures that the two-body resonances remain 
unaffected within the three-body system after this necessary fine-tuning. The 
effective total potential entering is then given by Eq.(4). 

In the outer part of Fig. lb the curves close to the horizontal zero axis show 
three gaussian three-body forces with different ranges and strengths adjusted 
to reproduce the measured two-neutron separation energy of 0.94 MeV for 
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Table 8 

The 2"'^ and 3'"'^ columns give the experimental and computed bound state (l/2~) 
and decay energies (3/2", 5/2", 7/2~, 9/2") in I'^Ne (in MeV). The 4*'' to 6*'* 
columns are the strengths of the gaussian three-body forces that for a range of 
3.0 fm, 3.5 fm, and 4.0 fm, give rise to computed energies in agreement with the 
experimental values. The 7^^ column is the expectation value of the three-body force 
with range of 4.0 fm. The last column is the contribution to the norm of the first 
three terms in the expansion (1) using the three-body potential of range 4.0 fm. 



J"" 


E 


^comp 


3.0 fm 


3.5 fm 


4.0 fm 


(Fab) (MeV) 


\i=l,2,3 (%) 


1/2- 


-0.94 


-0.79 


-1.5 


-0.9 


-0.6 


-0.2 


88.5, 11.1, 0.4 


3/2- 


0.34 


0.63 


-3.4 


-2.0 


-1.4 


-0.3 


90.7, 8.9, 0.2 


5/2- 


0.82 


0.91 


-1.1 


-0.7 


-0.4 


-0.1 


77.2, 16.9, 5.6 


7/2- 


2.05 


2.24 


-1.8 


-1.1 


-0.8 


-0.2 


97.5, 2.3, 0.2 


9/2- 


2.60 


2.70 


-0.1 


-0.1 


-0.1 


-0.1 


91.9, 4.2, 3.8 



the ground state of ^"Ne. The corresponding effective total potentials given by 
Eq.(4) are also shown in Fig. lb, but they can not be distinguished from each 
other within the thickness of the curves. This is because the three-body force 
is very weak compared to the depth of the full potential and it is furthermore 
largest for small p-values, where the total potential is highly repulsive. It 
is then clear that the main structure of the system can not be significantly 
modified by the choice of one or another of such three-body interactions. In the 
inner part of Fig. lb we show the computed effective potentials (4) for the five 
computed states of ^^Ne. The depth of the potential decreases with increasing 
total angular momentum. Higher J-values correspond then to higher energies 
of the ^■^Ne state. 



3.3 Spectrum of Ne 



The angular eigenvalues are computed from the effective potentials and used 
in the radial equations (3). Bound states are obtained as solutions falling 
off exponentially at large distances. Resonance eigenfunctions are found in 
complete analogy as exponentially falling solutions to the similar equations 
obtained after complex rotation of the hyperradius. 

Only the 1/2- state is bound in ^''Ne with a two-proton separation energy of 
—944 keV. The excitation energies of the four lowest resonance states are [13] 
1288±8 keV, 1764±12 keV, 2997±11 keV, and 3548±20 keV for the 3/2-, 
5/2-, 7/2-, and 9/2- states, respectively. These excitation energies corre- 
spond to the decay energies (energies above threshold) shown in the second 
column of table 8. The computed states (third column) are all slightly un- 
derbound compared to the experiment, but the ordering in the spectrum is 
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correct. Additional small attractive three-body potentials are needed to re- 
produce the experimental energies. For gaussians we give in table 8 both the 
necessary strengths for different ranges and the expectation values {Vsb{p)) 
for the corresponding J'^ solutions. The efficiency of the adiabatic expansion 
in Eq.(l) is seen in the last column of the table. Already the lowest adiabatic 
potential accounts for more than 75% of the probability and in most cases for 
more than 90%. The contributions from the different three-body potentials 
can hardly be distinguished. 

For the ground state of ^^Ne, we give in the last column of table 3 the per- 
centage of the total norm provided by each of the components defined by the 
sequence of quantum numbers in the previous columns. These numbers are 
computed with the three-body potential with range 4 fm, and they do not 
change significantly when another choice is made. When writing the three- 
body wave function in the first Jacobi set (upper part of the table) the main 
contribution of 89% is from s- waves, and less than 4% comes from rf-waves. If 
we use any of the other two identical Jacobi sets (lower part of the table) the 
d-waves contribute with roughly 50% while 45% is from the s-waves and the 
remaining few percents come from p-waves. 

It is remarkable that the interference between s and d waves explicitly included 
in the calculation plays a negligible role and falls below the 1% limit included 
in the table. The reason is that rotation of the components from one Jacobi 
set into another preserves the value of the total orbital angular momentum L. 
Then sd terms {L—2) in the second and third Jacobi sets must correspond to 
the almost non contributing {ix=£y = 1- L=2, 3=3/2} component in the first 
Jacobi set, since L=2 otherwise only arises from ix=0 or 2 where the anti- 
symmetry dictates zero spin of the two protons. Then L=2 and the resulting 
total spin of 3=1/ 2 can not couple to the total angular momentum 1/2 of the 
ground state. Thus this type of interferences is essentially excluded. 

The unbound excited states of ^^Ne are computed by application of the com- 
plex scaling method. The resonances obtained for ^^Ne are extremely narrow 
with widths much smaller than the accuracy of our calculations. Actually in 
[14] an experimental lower limit of 26 ps on the lifetime of the two-proton 
decay of the 3/2~ state is given. This lower limit for the lifetime amounts 
to an upper limit of 2.5 ■ 10~^^ MeV for the width of the resonance. Thus, 
application of the complex scaling method allows the use of very small scaling 
angles. Typically complex scaling angles of 6' = 10~^ are able to find the ^''Ne 
resonances. For these scaling angles the complex scaled A's can hardly be dis- 
tinguished from the non-rotated functions in Fig. la. The imaginary parts are 
very small and would appear on the zero line if plotted on the figure. 

In tables 4— 7 we show the largest components included in the calculation of 
the 3/2-, 5/2-, 7/2", and 9/2" states, respectively. For the 3/2" and 5/2" 
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states (as for the l/2~) the s, p, and d-waves are enough to get an adequate 
description of the system. Inclusion of f > 2 waves and use of extrapolated 
potentials does not significantly modify the results. The s-wave interaction 
between proton and core enters only in combinations with the d-wave. The 
component {£x—0, ^y—^-, L—0 } is simply not allowed in the 5/2~ state (the 
three intrinsic spins of 1/2 can not couple to 5/2). For the 3/2~ state rotation 
of the component {ix=iy=L=0, Sx=l, 3=3/2} from the second and third 
Jacobi sets into the first one contributes very little through the coupling of 
ix—i and iy=l due to antisymmetry and angular momentum conservation. 

Contrary to the ground state, the srf-interferences contribute substantially for 
both the 3/2^ and 5/2~ wave functions. For the 3/2^ state they give more 
than 70% of the probability in the first Jacobi set, and about 80% in the other 
two, see table 4. For the 5/2~ state these probabilities are roughly 73% in the 
first Jacobi set and 77% in the second and third, see table 5. In both states 
the j9- waves have a significant contribution only in the first Jacobi set, while 
in the other two sets all the possible {£x— ^j/=l} components give less than 
3% of the probability. 

For J^=7/2~ and 9/2~ (tables 6 and 7) s-waves do not contribute if only s, p, 

and (i-waves are included. For J'^=9/2~ the S(i-interferences are not allowed. 
For the 7/2^ state in principle the sd-interference could contribute only in the 
second and third Jacobi sets (L=2 and S=3/2 can couple to 7/2). However, 
rotation into the first Jacobi set preserves both the L=2 and 3=3/ 2 values 
implying Sx=l and Ix—i due to antisymmetry. These components give a con- 
tribution to the norm smaller than 1%. For these two resonance states only 
proton-core (i-waves for both Ix and ly are significant and appearing in the 
lower part of tables 6 and 7. 

When higher waves are included they do not contribute above the 1% level in 
the second and third Jacobi sets. The reason is simply that it is too expen- 
sive to increase the centrifugal barrier. However, the components in the first 
Jacobi set follow by rotation. In particular the large components with total 
orbital angular momentum L=3 and 3=3/2 must correspond to equally large 
contributions of odd partial waves in the first Jacobi set, which necessarily has 
£x > 3. For 9/2^ this is in fact the dominating component. The same happens 
with the L=4 and S'=l/2 components, that after rotation into the first Jacobi 
set give important contributions corresponding to s^f-interferences. 



3.4 Spatial distribution 

In Fig. 2 we show the square of the radial wave functions in the expansion (1). 
This allows a direct comparison of the contribution of each term to the total 
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Fig. 2. Outer part: Square of the radial wave functions corresponding to the two 
lowest effective potentials A„(p) for the l/2~ (sohd curves), 3/2~ (short-dashed 
curves), and 5/2~ (long-dashed curves) states of ^''Ne. Inner part: The same as the 
outer part but for the lowest effective potential A„(p) of the 7/2~ (solid curve) and 
9/2"" (dashed curve) states. 




Fig. 3. Contour diagram for the probability distribution of the l/2~ ground state of 
^^Ne. The square of the three-body wave function is integrated over the directions 

of the two Jacobi coordinates. 

wave function. In the outer part of the figure we show the first two terms in 
the expansion for the states with 7^^=1/2^, 3/2^, and 5/2~. For 7^^=7/2" and 
9/2^, shown in the inner part, only the contribution of the dominating term 
in Eq.(l) is visible. 

The complex scaling method provides bound states and resonances simulta- 
neously as solutions of the complex rotated equation (3) that asymptotically 
decrease exponentially. The radial solutions in Fig. 2 show this asymptotic 
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Table 9 

For the different computed states of -"^^Ne, root mean square distances {rpp}^^"^, 

('^c.pp)^^^' (^cp)^^^' (^Pjcp)"*^^^! (P^)^^^) aiid (r^)^/'^, where]? and c denote a halo proton 
and the core, respectively. All the distances are given in fm. The results do not 
depend on the three-body interaction chosen from table 8. 







/ 2 U/2 
\ c,pp/ 




/ 2 u/2 
\ p,cp/ 


(p2)l/2 


(r2)i/2 


1/2- 


4.5 


3.2 


3.9 


3.9 


5.3 


2.8 


3/2- 


5.3 


3.4 


4.3 


4.4 


5.9 


2.9 


5/2" 


5.5 


3.5 


4.3 


4.4 


5.9 


2.9 


7/2- 


5.6 


3.5 


4.3 


4.5 


6.0 


3.0 


9/2" 


5.7 


3.6 


4.4 


4.5 


6.1 


3.0 



behaviour already at distances of around 9 fm. 

For the excited states of ^^Ne, computed by the complex scaling method, the 
imaginary parts of the radial functions are very small and only the real parts 
are visible in the figure. This is a direct consequence of the extremely small 
widths of these resonances which allow a very small complex scaling angle 
resulting in almost real wave functions. As for any resonance, the true ra- 
dial wave function is non-square integrable, diverging at large distances as 
exp {\k\p sin [3 ji), where f3ji is the argument of the resonance. Due to the ex- 
tremely small values of (3r for the present cases of ^^Ne this divergence has no 
effect for p-values up to several hundred fm. 

Therefore, up to distances around 40-50 fm we can treat the resonances as 
bound states, and compute observables like (p^), (x^) or (y^), that give infor- 
mation about the spatial distribution of the three-body system. For each state, 
we give in table 9 such characteristic root mean square distances between the 
particles. The root mean square radius of the whole three-body system (r^)^/^ 
in the 1/2" ground state is 2.8 fm consistent with the experimental value of 
2.75±0.07 fm given in [32]. 

The structures derived from these mean values are crudely speaking triangles 
where the distance between the two protons is slightly larger than the distance 
between proton and core. The dimension increases with angular momentum. 
A more detailed view is obtained in the contour diagram in Fig. 3. We plot the 
square of the three-body wave function multiplied by the phase space factors 
and integrated over the directions of the two Jacobi coordinates. We observe 
that a prominent peak is found for distances between the two protons of about 
2 fm and the two-proton center of mass distance to the core of about 3 fm. 
Another smaller peak is found for corresponding distances of about 5 fm and 
1 fm, respectively. These two shapes can be described as two protons either 
on the same side of the core (at some distance apart) or at almost opposite 
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sides of the core. The structures of the excited states are very similar apart 
from the dimension increasing with J. 

The distribution overlaps to some extent the core density with the potential 
danger of violating the basic few-body assumptions. However, an assessment of 
model validity is only possible by comparing many predictions with measured 
observables. 

The size of a three-body system is for quantum halos related to the binding 
energy. The repulsive Coulomb potential may provide a very small binding 
energy but at the same time confining the structure spatially. The character- 
istic halo features arc revealed by the hyperradial extension measured in units 
of an appropriate length scale po, designed in [8] to measure the content of 
non-classical probability. In our calculation we have po ~ 5 fm, which implies 
that the dimensionless measure of the binding energy mBpl/h^ is around 0.6 
and the mean square hypcrradius in units of Po is about 1.2. This wave func- 
tion is therefore to a large extent in the classically allowed region and the 
properties cannot be obtained by scaling relations independent of the details 
of the interactions. The ^^Ne system is not a quantum halo as defined in [8] 
although the structure still can be described as a three-body system. 



3.5 Widths of the resonances 

The positive energies of the excited states imply that they can decay into three 
separate particles and also to lower-lying states by electromagnetic transitions. 
However, the widths of the resonances are very small and beyond the accu- 
racy of the otherwise very efficient complex scaling method. Fortunately the 
hyperspherical adiabatic method provides generalized radial potentials which 
are responsible for the three-body bound state structure as well as the asymp- 
totic behavior of the continuum wave functions. The decay widths into sepa- 
rate free particles can then be estimated in the WKB approximation. This is 
especially tempting since about 90% of the wave function is described by the 
lowest adiabatic potential and since the widths are small and the effective bar- 
rier is smooth. In [31] the accuracy of the WKB approximation was estimated 
for a system of three charged particles. The WKB turned out to be accurate 
to within a factor of five including the uncertainty from the preexponential 
factor, the "knocking rate" . 

This use of the WKB approximation, where the lowest adiabatic potential pro- 
vides the tunneling barrier, is a rather subtle application. There is a striking 
resemblance with computation of «-decay widths from tunneling through cen- 
trifugal and Coulomb barriers. However, a-decay is assumed to be a two-body 
process where the a-particle is formed with some probability and then tunnels 
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Fig. 4. Barrier of the most contributing effective potential for the 7^^=3/2", 5/2~, 
7/2~, and 9/2~ resonances in -"^^Ne. The crosses show the turning points and the 
arrows show the energies corresponding to each of the resonances. 

through the two-body barrier. The two-proton decay is a three-body problem 

and it is only due to the efficiency of the hyperspheric expansion that one 
suitable generalized potential becomes available. The tunneling probability 
through the corresponding barrier is, in the one A-approximation, then pre- 
cisely describing the resonance width for decay into the corresponding channel 
of three particles in the final state. On this level there is a complete equiv- 
alence with the a-decay computation. The WKB expression is then a rather 
good approximation to the full three-body coupled channel problem under 
the two conditions that the lowest adiabatic hyperspheric potential and the 
WKB-approximation both provide sufficient accuracy. 

We therefore have to compute the transmission coefficient through the barrier 
of the lowest effective potential V^q (see Eq.(4)) constructed by using the 
dominating angular eigenvalue, see Fig. lb. If crossings of potentials occur 
along the way through the barrier we follow the path determined by minimum 
change of the wave function. This corresponds to a path along the most smooth 
adiabatic potential in full agreement with the basic assumptions of the WKB 
approximation. Thus the WKB transmission coefficient is estimated by 



f^r [2,771 
-2 J [-^{yeSip) - Er) 



2m - -1'/^ . 

dp } , (7) 



where Er is the energy of the resonance, and pi and po are the inner and outer 
classical turning points defining the distance through the barrier. In Fig.4 we 
show the barriers of the resulting effective potential for the four resonances 
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Table 10 

Turning points, frequency, transmission coefficients, average lifetimes, and reso- 
nance widths for the 3/2~, 5/2", 7/2~, and 9/2~ resonances in -"^^Ne. pi and po are 
given in fm, / in MeV/^, ta in s, and V in MeV. 





Pi 


Po 


/ 




T 






ta 






r 




3/2- 


6.4 


46.7 


3.90 


9.2 


10- 


13 


1.8 


■ 10- 


-10 


3.6 


10- 


12 


5/2- 


6.4 


31.3 


3.82 


3.3 


10- 


11 


3.6 


• 10- 


-12 


1.3 


10- 


10 


7/2- 


5.5 


15.4 


1.54 


3.7 


• IQ- 


-3 


1.2 


• 10- 


-19 


5.6 


• IQ- 


-3 


9/2- 


5.5 


13.7 


1.59 


3.3 


• 10" 


-3 


1.2 


• 10- 


-19 


5.3 


• 10" 


-3 



in -"^^Ne with J'^ = 3/2-, 5/2-, 7/2- and 9/2-. The outer turning points of 
the effective potentials are located at distances of around 15 fm for the 7/2- 
and 9/2" states, above 30 fm for the 5/2" resonance, and above 45 fm for 
the 3/2- resonance. These potentials, and therefore the angular eigenvalues, 
must then be computed accurately below these values. Inaccuracies typically 
produce too high eigenvalues and therefore too large barriers with too small 
transmission coefficients or equivalently too small widths. 

Once the transmission coefficient is computed, the decay constant can be ob- 
tained as V/h—fT, where the frequency / is given in terms of the radial 
extension pi of the attractive pocket and the velocity v of the particle moving 
in it, i.e. f=v/pi. We estimate the velocity v by equating the kinetic energy 
{niv'^/2) and the maximum potential energy E'/j+l V^^g*""* |, where |V|fg*"^| is 
the maximum depth of the effective potential. This estimate of the knocking 
rate can easily be off by a factor of two. Furthermore this computation implic- 
itly assumes that the preformation factor is unity, i.e. that the decaying state 
really is a three-body state consisting of two protons outside the ^^O core. 

Finally, the average hfetime ta or the width F of the resonances are found 
by inverting the decay constant or multiplying by h, respectively. In table 10 
we summarize the results for the different resonances. The computed width 
of the 3/2- resonance is consistent with the upper limit of 2.5 ■ 10"^^ MeV 
quoted in [14] . This limit is in fact essentially theoretical as obtained in a shell 
model calculation for the magnetic dipole deexcitation probability. The lack 
of protons in the experiment then leads to the conclusion that the width must 
exceed the computed electromagnetic value. 

Similarly a lower limit for the two-proton decay of the 5/2" resonance is de- 
duced to be about 3 • lO-^'^ MeV [14] which is about three times larger than 
our crude calculation shown in table 10. The two higher-lying resonances have 
much larger widths of about 5 keV due to the higher energies and the rather 
narrow barriers. This behavior can not be predicted from angular momen- 
tum and energy systematics. The three-body model is indispensable, not only 
because the final state asymptotically consists of three particles but more in- 
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terestingly because the short and intermediate distances determine the widths. 



4 Accuracy of model results 

The calculated results can not be more rehable than the model. The basic 
assumptions are that the active few-body degrees of freedom are responsible 
for all the computed properties. The intrinsic constituent particle degrees of 
freedom are only implicitly treated through the properties of the two-body in- 
teractions. Polarization effects are included but not explicit contributions from 
intrinsic structure. Therefore cluster models are most suited in descriptions 
of halo systems while individual applications beyond halos still can be rather 
accurate if the different degrees of freedom are sufficiently weakly coupled. 

Ultimately only practical tests can decide on the issue of model rehability. 
However, a necessary prerequisite is a series of accuracy tests of the model 
itself. This means that the numerical methods must reach convergence to the 
requested level of accuracy. This is not trivial in the present type of cluster 
computations. Another source of uncertainty is the parametrization of the 
interactions, e.g. spin-dependence, radial form factors and the numerical values 
of the parameters. We shall in this section report on these tests. 

4.1 Convergence and basis expansion 

The convergence of the calculations must be reached at three different lev- 
els. First, in the hyperspheric adiabatic expansion of Eq.(l) the number of 
included potentials must be sufficient. This was shown in the last column of 
table 8. Typically three terms are enough to provide contributions to the wave 
functions larger than 99% of the probability. Sometimes one more adiabatic 
potential has to be included for additional confidence. 

Second, the expansion of the angular eigenfunctions 0^*^ in Eq.(2) in terms of 
hyperspherical harmonics has to contain all the components giving a relevant 
contribution to the wave functions. These are the components shown in ta- 
bles 3 to 7. Together with them, additional components giving contributions 
smaller than 1% have also been included in the calculations. 

Third, once the components have been chosen, convergence must be secured 
in the expansion in hyperspherical harmonics by requiring inclusion of a suf- 
ficiently large number of them, i.e. y^^^ i{ai,flxi,^yi)- For each component 
we include all these basis functions with hyperspherical quantum number less 
than a maximum number K^ax which may depend on the quantum numbers 
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Fig. 5. Comparison of the computed A functions for different values of K^ax for the 
i^Ne states J"" = 1/2" (left) and J'^ = S/2~ (right). The solid curves are the results 
with K^ax from tables 3 and 4. 

of the components. We should emphasize that the basis sizes we discuss are 
for each Faddeev component, i.e. related to each Jacobi coordinate system. If 
only one Jacobi set is used a much larger K^ax is necessary for convergence. 

In the calculations this hypermoment Kmaxi given in the sixth column of tables 
3-7, is by far large enough to reach convergence in the region of contributing 
p- values. This is illustrated in Fig. 5 where we show the three first angular 
eigenvalues, the A-functions, for the l/2~ and 3/2~ states in -"^^Ne for three 
different values of Kmax- For the ground state we observe that at a distance 
of around 7 fm the lowest A with Kmax=^^ begin to differ from the converged 
result. From Fig. 2 we see that this significantly can modify the tail of the 
l/2~ radial wave function. In fact, using the three-body interactions given in 
table 8 the two-neutron separation energy of the l/2~ state with Kmax=^^ is 
only —0.58 MeV, almost 400 keV less bound than in the full calculation. For 
Kmax='^^ the discrepancy from the converged lowest A for J'^=l/2~ begins at 
distances larger than 10 fm, and the effect on the radial wave function is not 
important. When K^ax—'^^ the two-neutron separation energy is —0.93 MeV, 
pretty close to the one we obtained. Therefore for the ground state at least a 
value of Kmax='^^ is needed. 

For the first excited state, 3/2~, the discrepancies between the different calcu- 
lations begin at lower distances, see Fig.5. When Kmax=^0 the lowest A begins 
to differ from the other calculation already at p=6 fm. At this value of p the 
radial wave function in Fig. 2 is close to its maximum, and therefore the results 
obtained with Kmax=^0 are not reliable. Actually, the 3/2~ resonance energy 
obtained in this case is above 700 keV, more than a factor of two higher than 
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in the full calculation. For K„iax='20 the computed 3/2~ resonance energy is 
still 15% higher than the 340 keV obtained in the full calculation. The same 
kind of convergence is seen for higher values of J. 

We can then conclude that for the excited states of ^''Ne values of K^ax higher 
than 20, typically around 30, are desirable. As shown in tables 4-7 we used 
in our calculations Kmax values significantly larger than 30 for some of the 
components. This is because the estimates of the widths of the resonances 
require accurate effective potentials for much larger distances than 10 fm, e.g. 
up to about 50 fm for the 3/2~ resonance. This accuracy is achieved with the 
Kjnax values given in the tables. Thus reliable width calculations in a three- 
body model require a rather large individual basis for each of the Faddeev 
components related to the different Jacobi sets. 

4-2 The shape of the radial potential 

In this section we investigate the dependence of the results on the radial shape 
of the nuclear two-body interactions. The results discussed so far have been 
obtained with the proton-proton interaction (G) from [26,29] and the Woods- 
Saxon (WS) proton-core interaction in Eq.(5). 

Let us consider now also the gaussian (G) proton-core interaction given in 
Eq.(6) and table 1 and the sophisticated Argonne (A) proton-proton interac- 
tion described in [33]. This is a nonrelativistic potential that reproduces 
pp and np scattering data for energies from to 350 MeV, nn low-energy 
scattering data, and the deuteron properties. Different combinations of these 
potentials permit comparison and reliability tests of the results. The calcula- 
tions are denoted G+WS, A+WS, G+G, where the first and second labels 
refer to the proton-proton and proton-core interactions, respectively. 

We have repeated all the three-body calculations for the G-l-G and A-I-WS 
cases. The decisive A-functions entering in the radial equations (3) are com- 
pared in Fig.6 for the 1/2^, 3/2^, and 5/2~ states in ^''Ne. At large distances 
the A-functions are basically identical in the three calculations, since the low 
energy properties of the two-body potentials are the same. At short distances 
some differences are visible. The G+G calculations produce a slightly less deep 
pocket in the potentials than in the other two cases. This means that the phe- 
nomenological three-body potentials needed to fit the measured energies of 
the ^^Ne states are also slightly different. 

In total, after fitting the resonance energies in ^''Ne by use of gaussian three- 
body forces, the results are essentially the same as the G+WS results previ- 
ously discussed. The radial wave functions are hardly distinguishable from the 
ones shown in Fig.2. In Fig. 7 we compare the weights of the components for 
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Fig. 6. Comparison between the deepest A-functions for the 1/2^, 3/2~, and 5/2~ 
states in ^^Ne when different proton-proton and proton-core interactions are used 
(G=Gaussian, WS= Woods-Saxon, A=Argonne). In the label the first symbol refers 
to the proton-proton potential and the second to the proton-core. 

the three calculations and the five computed states in ^''Ne. The components 
are ordered as in tables 3 to 7. The left part of the figure corresponds to the 
components when the full three-body wave function is written in the first Ja- 
cobi set, and the right part to the ones in the second and third Jacobi sets. 
The three computations give very similar results. For each state the dominant 
components are the same in all the three cases, and the overall three-body 
structure is essentially unchanged. 

We can then conclude that details like the radial shapes of the two-body 
interactions are not changing the properties of the three-body system within 
the accuracy provided by the inherent accuracy of the model. 



5 Previous works 

Different theoretical works concerning the structure of ^^Ne can be found in 
the literature. Most of them [15,17,18] focus only on the ground state structure 
and in particular on the weight distribution of the s and (i-waves. Only [16] 
and [19] report on more detailed calculations of the ^^Ne spectrum. 

In [16] the three-cluster generator coordinate model is used to obtain ^''Ne. The 
proton-proton interaction is adjusted to reproduce the experimental binding 
energy. Although no numbers are quoted, the s-wave contribution is expected 
to dominate. The same conclusion is reached in [17] where they analyze the ex- 
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Fig. 7. Comparison of the weights of the components for the ^^Ne states. The left 
part refers to the first Jacobi set and the right to the second and third Jacobi sets. 
The components are ordered as in tables 3-7. For each component the three bars, 
from left to right, refer to the G+WS, A+WS, and G+G computations, respectively. 

perimental parity dependence of the Coulomb shifts of the A=17 isoquartets. 
They used harmonic oscillator wave functions with a state dependent radius 
and the Coulomb potential of a uniformly charged sphere. The two protons of 
the ground state of ^''Ne then seem to occupy the si/2-orbits. 

In [18] the ^^-configuration is for the first time found to dominate in the 
ground state of ^''Ne. This is concluded by investigating the (3^ decay to the 
first excited state of ^^F. The decay is found to be roughly a factor of two 
larger than expected from the nuclear matrix elements that reproduce the (3~ 
decay into ^'''N. In [15] the domination of the ^^-configuration is confirmed 
by investigating the difference in Coulomb energy for s and ti-orbitals. They 
use Woods-Saxon potentials and an uniformly charged sphere. A coupling of 
either two Si/2 or two (i5/2-protons to the A=15 core immediately reveals that 
the (i^-component dominates. To reproduce the measured ^''Ne two-proton 
separation energy around 78% of the wave function must be of d^-character. 
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Table 11 

Computed s, p, and d wave contents in the ^^N (upper part) and ^''Ne (lower part) 
wave functions. The results in parenthesis are the ones obtained in ref.[19]. 



J"" 


s p d 


3/2- 
5/2- 


38 9 f39 8) 6 7 (4 5) 54 4 (55 6) 

32.8 (36.8) 2.2 (3.0) 65.0 (60.0) 

32.9 (34.5) 4.2 (3.2) 62.9 (61.6) 


1/2- 
3/2- 
5/2- 


45.4 (48.1) 5.6 (4.0) 49.0 (47.8) 
39.2 (38.1) 2.1 (2.8) 58.7 (58.4) 
37.4 (36.2) 3.3 (3.0) 59.3 (60.1) 



A detailed calculation of -"^"^Ne and ^'^N has been reported recently in [19]. A 
three-body model is employed and Woods-Saxon potentials between the three 
pairs of particles are used. The weight of the different partial waves for the 
ground state contains around 40% and 48% s-wave in the ^^N and ^^Ne wave 
functions, respectively. These values decrease for the first and second excited 
states, where the s^-wave content varies from case to case between 34.5% 
and 38.1%. In table 11 we compare our results to the ones given in [19]. For 
completeness we also include our results for the ^^N case. The agreement is 
rather good in all the cases. However, apparently in [19] interferences between 
s and d waves are not considered, while we have found that for the excited 
states they are crucial. Furthermore, it is not easy to understand how two 
protons in the si/2-waves can give rise to the excited 5/2- state, provided the 
^^O core still has spin of 1/2. It is also important to remember that by £-wave 
content we refer to the weight of the components in tables 3—7 containing 
that value of i. In this way, our ix=('y=0 components are not fully equivalent 
to the s^-content in the shell model sense. 

The spectrum of ^'''Ne is previously investigated only in [16] and [19]. In [16] 
the nucleon-nucleon interaction is adjusted to reproduce the ^^Ne ground state 
binding energy, but the same interaction is not able to describe the ^^N and 
^^F spectra. Even the level order is wrong and the 3" level comes out as 
the ground state in both cases. For the three-body systems the ground state 
angular momentum is correct in [16], but the 3/2- and 5/2- levels are reversed 
compared to the experimental data. The lack of consistency between the two- 
body interaction needed to reproduce the three-body ground state and the one 
reproducing the two-body spectra of the two-body subsystems is certainly a 
weak point in this calculation. 

The same deficiency is present in [19], where the three-body structure in ^'''Ne 
obtained with the two-body interactions that reproduce the ^^F spectrum 
also reverse the positions of the 3/2" and 5/2--states. The 5/2--level is even 
bound by —0.41 MeV, while the 3/2--state has an energy above threshold 
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of 0.89 MeV. A three-body force adjusted for each of the states is needed to 
restore the level order and the positions as found in experiments. It is then 
clear that the three-body force is forced to play a too important role, since 
the 5/2~ energy is moving up from —0.41 MeV to the experimental 0.82 MeV 
and simultaneously the 3/2 ~ is moved down in energy. 

In the present work these problems are not encountered. When only the two- 
body forces describing the proton-proton and ^^F properties are used, the 
energies of the ^^Ne states are the ones in the third column of table 8. They 
follow the experimental ordering, and use of a small effective three-body force 
is enough to fit the experimental energies. 

The fact that in [19] the same inversion as in [16] appears can be related to 
the inappropriate spin dependent operators used in the two-body nucleon- 
core interaction. Contrary to the present work the interaction in [19] is mixing 
the and the (i5/2-states used as an essential part of the valence space for 
the outer nucleons. The low-lying resonances in ^^N and ^^F do not have the 
desirable pure (i5/2-character [28]. 

Finally the partial two-proton decay widths of the three-body resonances are 
estimated in the present work, see table 10. In [19] the widths of the 3/2~ 
and 5/2~-states are found to be 4.1 ■ 10^^^ MeV and 1.2 ■ 10^"*^^ MeV, respec- 
tively smaller by seven and one order of magnitude compared to our results. 
This computed 3/2^ value is in apparent agreement with the experimental 
limit. However, their computed width of the 5/2~ resonance is already at 
least one order of magnitude too small. The eight orders of magnitude differ- 
ence between these two computed values seem to be in clear disagreement with 
measurements which indicates that the lowest and highest of these resonances 
have widths slightly below and above the electromagnetic decay widths. 

The origin of these discrepancies can not be decided without full repetition 
of all computations. Two effects are perhaps important in this connection. 
The first is that the basis size has to be rather large to provide a sufficiently 
accurate description all the way up to the relatively large distances where 
important contributions to the widths are determined. Too small a basis would 
usually lead to too high potential energies and therefore too large barriers and 
too small widths. The maximum hypermoment must be very large when only 
hyperharmonics in one Jacobi set are used. 

The second effect is that the hypersherical adiabatic expansion in principle 
includes all decay channels, i.e. also sequential decay which often is the domi- 
nating mode. For the 3/2~-state this is apparently not allowed due to energy 
conservation. However, a kind of virtual decay, where the tail of the two-body 
resonance is exploited, would enhance the decay probability. This effect is 
included in our formulation since the two-body resonance configuration influ- 
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ences the lower adiabatic potentials to fairly large hyperradii. This effect is 
difficult to catch with the direct hyperharmonic expansion method. 



6 Summary and conclusions 

The structure of the lightest Borromean nucleus ^^Ne is investigated in details 
in a three-body model where two protons surround the ^^O core. We use the 
hyperspheric adiabatic expansion method with two-body interactions adjusted 
to reproduce the properties of the two-body subsystems. The spin-dependent 
form of the proton-core interaction must be consistent with the mean-field ap- 
proximation. Otherwise the assumed decoupling of core and valence motion is 
violated due to the Pauli principle and the angular momentum conservation. 
We use two different radial shapes for the proton-core potential both consis- 
tent with the measured size of ^^O. The Coulomb potential is derived from a 
gaussian charge distribution with the measured root mean square radius. 

The three-body ground state and four measured excited states of ^''Ne are 
computed. The structure for the ground state is found to be about 45% of 

and 50% of (P {i^—^y—'^) proton-core components. For the two low- 
est excited states the sd relative proton-core states are dominating while only 
roughly 18% is of rf^-character. For the two highest resonances, (i^-components 
dominate completely in the second and third Jacobi sets, while contributions 
from higher partial waves {sg and ^/-components) can be important in the 
first Jacobi set. Properties of these states are not computed in previous works. 

The spatial distribution is slightly more extended than for ordinary nuclei and 
with a clear tendency to show two components, i.e. one where the protons and 
the core are distributed in a triangle and one where the protons try to be on 
opposite sides of the core. The distribution is not in the classically forbidden 
region and the scaling relations specific for quantum halos are not obeyed. 
The criteria for quantum halos are not fulfilled. Still the three-body model 
seems to give a rather good description of the structure of these five lowest- 
lying states. The two-proton decay widths are computed by use of the WKB 
approximation applied on the adiabatic potentials. The computed results are 
within the experimentally determined limits. The width estimates employ a 
recently established but non-trivial formulation which relies on the efficiency 
of both the hyperspheric adiabatic method and the separation at small and 
intermediate distances of the intrinsic and the necessary asymptotic three- 
body degrees of freedom. 

Accuracy of these computations can be divided into three parts, i.e. (i) model 
assumptions related to its applicability, (ii) input determination, uniqueness 
and reliability, (iii) numerical accuracy including convergence of expansions. 
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First the consistency of the results supports the crcdibihty of the modeL Sec- 
ond we carefuhy evahiated the accuracy both by testing for convergence and 
by using different two-body interactions. 

In conclusion, the three-body model describes efficiently the cluster structure 
of ^''Ne. The hyperspheric adiabatic expansion in the same framework provides 
consistently also the simple semiclassical estimates of resonance decays into 
three-body final states. 
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